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Abstract. Generalized Parton Distributions of the pion are evaluated in chiral quark models with 
the help of double distributions. As a result the polynomiality conditions are automatically satisfied. 
In addition, positivity constraints, proper normalization and support, sum rules, and soft pion 
theorems are fulfilled. We obtain explicit expressions holding at the low-energy quark-model scale, 
which exhibit no factorization in the f-dependence. The crucial QCD evolution of the quark-model 
distributions is carried out up to experimental or lattice scales. The obtained results for the Parton 
Distribution Function and the Parton Distribution Amplitude describe the available experimental 
and lattice data, confirming that the quark-model scale is low, around 320 MeV. 
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Generalized Parton Distributions (GPD's) carry "tomographic" information on the 
partonic structure of hadrons (for reviews see e.g. [1, 2, 3, 4, 5, 6, 7, 8]). In this talk we 
present our recent calculation of the GPD's of the pion in the framework of chiral quark 
models [9], which extends the previous calculations of PDF's [10, 11, 12], PDA's [13, 
14], and GPD in the impact parameter space [15]. Recently, the Transition Distribution 
Amplitudes (TDA) [16, 17] have also been evaluated in the same framework [18]. 
Other quark- model calculations of GPD's and related quantities have been reported in 
Refs. [19, 20, 21, 22, 23, 24, 25, 23, 26, 27, 28, 29, 30]. 

Chiral quark models yield parton distributions at a given low energy scale Qq. The 
result for a quantity F is matched to QCD order by order in the twist expansion, n, hence 
^;!(^)|Modei = Fn{x, Oc{Q^))\Qciy- Then the functions F„ are evolved to higher scales Q. 
It turns out that in order to describe the available pion phenomenology the initial scale 
Qq in the considered quark models must be very low [10, 12, 13, 14]; matching the 
momentum fraction carried by the valence quark at = 4GeV^ to 47% [31, 32] yields 

eo = 313l?gMeV, (1) 
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FIGURE 1. The direct (a), crossed (b), and contact (c) Feynman diagrams for the quark-model evalua- 
tion of the GPD's of the pion. 



with Aqcd = 226 MeV and three flavors. At such a low scale a{QQ)/{2'Jl) = 0.34, which 
makes the evolution very fast for the scales close to the initial value. 

The kinematics of the process and the assignment of momenta (in the asymmetric 
notation) is displayed in Fig. 1, representing the large-A^^ quark-model evaluation of 



GPD's. We adopt the standard notation 
leading-twist GPD of the pion is defined as 
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where a and b are isospin indices for the pion, T is the isospin matrix equal 1 for 
the isoscalar and T3 for the isovector case, n is the null vector, and z is the light-cone 

coordinate. In the symmetric notation one introduces ^ = and X = jr^- The 
following sum rules hold on general grounds: 

jWj^^=\x,^,t) = 2Fy{t), jWxj^'=^{x,^j) = e2{t)~^^ei{t), (3) 

where Fv{t) is the electromagnetic form factor, while Oi{t) and 02(0 are the gravita- 
tional form factors of the pion. Finally, for X > the equality jr'^^'\X ,0,0) = q{X) 
relates the GPD's to the the pion's parton distribution function (PDF). The polynomiality 
conditions [1,2] and the positivity bound [33] are satisfied in our approach. 

We work for simplicity in the chiral limit, nin = 0. Two quark models are considered: 
the Spectral Quark Model (SQM) [34] and the NIL model. SQM implements the vector- 
meson dominance, predicting the form factors 
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FIGURE 2. The GPD's of the pion in SQM, B, = 1/3, f = 0.2,0, -0.2, -1, -10, -100 GeV^ from top 
to bottom (at x = 0.9). 



The explicit results for the full GPD's have been provided in Ref. [9]. Importantly, their 
form does not exhibit a factorized ^-dependence. A sample result for = 1 /3 and several 
values of t is shown in Fig. 2. For the NIL model the results are qualitatively the same. 
For the case of ? = the GPD's simplify to the well-know [20, 25] step-function results 

^=o(-^,C,0) = d[{\-x){x-0]-d[-x{x+l-Ol 
^=i(x,C,0) = e[{l-x){x+l-Q]. (5) 

Another simple case is in SQM for ^ = and any value of r [15] 

jE{x,0,t)= ^ ^ LL. (6) 

For the QCD evolution we use the leading-order ERBL-DGLAP equations with three 
flavors. In the left panel of Fig. 3 we confront the result for xq{x,Q) at the scale 
2 = 2 GeV with the data at this scale from the E615 Drell-Yan experiment [35]. We 
note agreement between the model and the data. In the right panel of Fig 3 we compare 
our results to the data from lattices [37]. We take the liberty of moving the scale, as its 
determination on the lattice is not very precise. As we see, the agreement is qualitatively 
good if one considers the uncertainties of the data, especially when the lower scale is 
used. 

PDA's have been intensely studied in the past in several contexts (see Ref. [38] for 
a brief review). At the quark model scale Qq the PDA of the pion [13], which can 
be related to the isovector GPD through the soft pion theorem [39] is ^{x;Qo) = 1 
[13]. The evolved PDA is shown in Fig. 4, where it is compared to the E791 di-jet 
measurement [40] and to lattice calculations [37]. Again, good agreement is observed. 

For the case of general kinematics, the explicit form of the LO QCD evolution 
equations for the GPD's can be found in [41, 42, 43, 44, 45, 46, 47]. In this paper 
we solve them with the numerical method developed in [45], based on the Chebyshev 
polynomial expansion. The results of the LO evolution from the SQM initial condition 
for = 1/3 at the scale Qo to subsequent values of Q are shown in Figs. 5. The evolution 
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FIGURE 3. Left: the quark model valence parton distribution (PDF) of the pion for a single quark 
(either u or d for ;r+) evolved to the scale of Q = 4 GeV (band). The width of the band indicates 
the uncertainty in the initial scale Qq. The data points come from the analysis of the E615 experiment 
[35]. The dashed line shows the reanalysis of the original data from Ref. [36]. Right; the quark-model 
prediction for PDF evolved to the scale 2 = 0.5 GeV (darker band) and 2 = 0.35 GeV (lighter band). The 
transverse-lattice data come from Ref. [37] and correspond to the scale ~ 0.5 GeV. The line shows the 
GRS parameterization at 2 = 0.5 GeV. 
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FIGURE 4. Left: the quark-model prediction for the pion distribution amplitude (PDA) evolved to the 
scale 2 = 2 GeV (band) and compared to the E791 di-jet measurement [40] after proper normalization 
of the data. The width of the band indicates the uncertainty in 2o- We also show the the asymptotic 
PDA, 0(x,oo) — 6x{l —x) (dashed line). Right: the same compared to the transverse lattice data [37], 
corresponding to the scale ^ 0.5 GeV. 



is fastest at low values of Q, where the coupling constant is large, and it immediately 
pulls down the end-point values to zero. Then, the strength gradually drifts from the 
DGLAP regions to the ERBL region. The approach towards the asymptotic form is very 
slow, with the tails in the DGLAP region present. The highest displayed in the figure 
is 10^ GeV^ and the asymptotic form is reached at "cosmologically" large values of Q, 
which are never achieved experimentally. The results for the NJL model are very similar 
to the case of SQIVI. 





FIGURE 5. Non-singlet (top), singlet (middle), and gluon (bottom) GPD's of the pion for ^ = 1 /3, 
evolved from the SQM condition up to = 0.1, 1, 10, 10^, . . . , 10** GeV^. Higher gives higher 
magnitude of the curves in the ERBL region. 



In conclusion, we remark that our calculation provides a link between the non- 
perturbative soft-energy physics in terms of matrix element of operators and the high- 
energy processes as deduced from perturbative QCD evolution. The overall agreement 
with the pionic data from experiments and lattices, available for the PDF and PDA, is 
very reasonable, supporting the presented methodology. 
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